Correlations in noisy Landau-Zener transitions 
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We analyze the influence of classical Gaussian noise on Landau-Zener transitions during a two- level 
crossing in a time-dependent regular external field. Transition probabilities and coherence factors 
become random due to the noise. We calculate their two-time correlation functions, which describe 
the response of this two-level system to a weak external pulse signal. The spectrum and intensity 
of the magnetic response are derived. Although fluctuations are of the same order of magnitude as 
averages, the results is obtained in an analytic form. 



PACS numbers: 75.50.Xx,05.40.Ca 

I. INTRODUCTION 

The Landau-Zener (LZ) theory (see also Q) plays 
an important role in many different physical problems 
ranging from chemistry, biology, theory of collisions to 
the tunneling of Bose-condensates, dynamics of glasses 
and spin reversal in magnetic wires and nanomagnets, as 
well as quantum computing. LZ theory treats a rather 
simple and general situation of two-level crossing. Near 
the point of (avoided) level crossing, adiabaticity is vio- 
lated and transitions occur. The LZ Hamiltonian can be 
represented by a 2 x 2 matrix acting in a two-dimensional 
space of vectors spanned by the basis vectors |f) and ||), 



Ht,z — 



£ T A 
A* E i 
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where Et 



ment A does not depend on time. Without loss of gener- 
ality one can assume A to be a positive constant. Landau 
and Zener have found the evolution matrix 



-tkvt and the transition matrix ele- 



Ulz 



a b 

-b* a* 



(2) 



which connects initial amplitudes in the diabatic basis Q 
with their final values. The entries 
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depend only on the dimensionless parameter 
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A 



(3) 



(4) 



which we call the LZ parameter. 

In many applications the LZ theory must be modified 
to take into account thermal noise or noise of a different 
nature. In a pioneering work Kayanuma |5| has solved 
such a problem on the basis of three simplifying assump- 
tions: i) the regular transition matrix element is equal 



to zero (i.e., A = 0), the noise is completely responsible 
for transitions; ii) the correlation function of noise has a 
simple exponential form; iii) the noise is fast. The first 
two limitations were lifted in the work by Pokrovsky and 
Sinitsyn [(|. Their analysis of time scales for different 
processes distinguishes the correlation time of noise r n , 
the accumulation time r acc = (i/T n )~ l during which the 
noise effectively produces transitions, and the LZ time 
tlz = A /(Hi/) during which the standard LZ transi- 
tion (without noise) proceeds. The noise is called fast if 
T n *C tlz *C T acc . In this case the separation of times al- 
lows to neglect the regular matrix element A at \t\ ^> tlz 
and to neglect the noise in the interval \t\ <C r acc . Since 
these intervals overlap, one can solve first the reduced 
problem with the noise producing transitions and then 
match this solution with the LZ solution. In this way the 
authors of Ref. |(| have found the average value of the 
density matrix. 

The two-level problem is as usual equivalent to the 
problem of spin i rotating in the time-dependent field 
containing regular and stochastic parts. In their second 
work [7| the same authors represented the density matrix 
in terms of the Bloch vector g, 



P(t) 



1 , 



[1 



(5) 



It is 



where <r is the triplet of Pauli matrices for spin i. 
a well known fact that the square of the Bloch vector is 
an integral of motion. In Ref. this fact was used to 
calculate the square fluctuation 



(mm = 



;(*)> 2 



However, the fluctuations of separate components, 
([5gi(t)] 2 ), (i — x,y,z) were not calculated, neither two- 
time correlators (gi(t)gj(t')) were found. 

Such fluctuations are observable, for example, in a sys- 
tem of magnetic molecules subject to the same mag- 
netic field (with identical regular and random contri- 
butions). The measurement of the z component of 
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the magnetic moment averaged over the system yields 
m z(t) = | Tr a z ■ p(t) — %g z {t) which depends on the 
particular noise history of the measurement. The repeti- 
tion of such a measurement for many independent noise 
histories then yields an average (m z (t)} with fluctuations 
([Sm z ] 2 ) = \{[5g z {t)]). Similarly, the twofold time de- 
pendence of the correlators (gi(t)gj(t')) determines spec- 
tral properties of the magnetization. Namely this situa- 
tion is realized in magnetic systems subject to an exter- 
nal regular and random magnetic field. In both cases we 
assume that the characteristic wave length of the electro- 
magnetic field is much larger than the linear size of the 
system. 

In two-level systems the random field can be realized 
by an environment, for example, by thermal phonons or 
by magnetic field of nuclei. In such a situation different 
two-level systems feel different random fields. Therefore, 
in the ensemble consisting of a large number N of the 
two-level systems the fluctuation will be suppressed pro- 
portional to 1/yN. However, the two-time correlation 
function describes the linear response of such a system 
to a small perturbation 6h(t) (identical for all members) 
according to a standard linear-response equation: 

S(ai) = f dt' {s t {t) Sj {t'))5h 3 (t'). (6) 

Such a perturbation can be realized as a pulse of electric 
magnetic or acoustic field. In this work we determine 
all these correlations under the same assumption of fast 
noise. 

The outline of the paper is as follows. In Sec. ITT1 we 
formulate the model for LZ transitions in the presence of 
noise we focus on. In Sec. IIIII we calculate in a general 
form the time evolution of averages and autocorrelations 
of the Bloch vector in the absence of the regular transi- 
tion matrix element A. These general findings are evalu- 
ated to obtain the averages and fluctuations of transition 
probabilities first for A = in Sec. II VI and then also for 
A > in Sec.[V] Implications for the spectral width and 
intensity of fluctuations for transitions in a gas of collid- 
ing atoms or molecules are discussed in Sec. I VII We close 
with concluding remarks in Sec. IV 111 Some calculational 
details are collected in an appendix. 



II. FUNDAMENTAL 

We consider the dynamics of a quantum spin S in 
the presence of a time dependent effective magnetic field 
B(i). Absorbing the Lande factor and the Bohr magne- 
ton into this field, the Hamiltonian reads 

H(t) = -B(t) • S . (7) 

The effective field B(t) = b(i) + r/(t) is composed of 
a regular part b and a noisy part r] with zero average, 
(ri(t)) = 0. In the LZ theory, an avoided level crossing is 



described by the regular field which may be chosen as 

b(t)=b z (t)e z +b x e x (8) 

with a z component approximately linear in time, 

b&)avt, u = b z {0), (9) 

and an approximately time independent perpendicular 
component b x = —(2/H)A. The noise is considered as 
Gaussian distributed and uncorrelated in the longitudinal 
and transverse directions with variances = x,y, z) 

(0) = <% A(*-0 • (io) 

The correlation functions fi naturally are even functions. 
A further planar symmetry f x = f y = f± is assumed. 

Under the action of the Hamiltonian @ , the dynamics 
of the density matrix is equivalent to the classical equa- 
tion of motion (Bloch equation) 

g(t) - -B(t) x g(t) . (11) 

For the analysis of the Bloch equation it is convenient to 
introduce 

9± = ~^(dx ± ig y ) (12) 

as well as T)±, b±, and B± by analogous definitions. In 
these terms the Bloch equation can be rewritten as 

g± = TiB z g± ± iB±g z , (13a) 
g z = i[B-g+ - B+g-} . (13b) 

From these equations the longitudinal field component 
B z (t) can be eliminated by going over to the interaction 
representation with respect to the Hamiltonian 

H (t) = -B z (t)S, . (14) 

Using the time evolution operator 

tr (t) =exp|i^ df#o(f)} (15) 
we rewrite the density matrix JSJ as 

p(t) = l^+u (m(t)-^ Q (t)} (i6) 

in terms of the transformed Bloch vector 

g±(t) = 3±(t)exp|±i^ dt'B z (t')^ . (17) 

The transformation does not affect the longitudinal com- 
ponent, g z (t) = g z (t). Using analogous transformations 
for 17, b and B, the Bloch equation simplifies to 

g± = ±iB±g z , (18a) 
~g z = i[B_g + - B + g_] . (18b) 
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In the limit of short noise correlation time r n , the cor- 
relations for the transformed noise read 

(fi + (t)fj + (t')) = <*7- (*>?-(*')> = 0, (19a) 
(r?+ (*)*?-(*')> = h (*,*') , (19b) 

with the correlator 

/j.(t,0 = / ± (i-i')exp|i J^dt"B z (t") 

« /x(t-t') e 4 ^* 2 "*' 2 )/ 2 . (19c) 



Since — i') decays for |t — i'| ^ r n one may neglect 
the contribution of rj z (t") to the time integral provided 



t 1 2 (^ 2 )«1 



(20) 



This condition is assumed in the following. It can be 
considered as a further limitation for the noise corre- 
lation time or as a limitation for the noise amplitude. 
However, it is not crucial and is introduced only for sim- 
plification. If the longitudinal noise is independent of the 
transverse one, its contribution to the averaged exponent 
in equation Ijl9(l can be reduced to the Debye- Waller fac- 
tor W(t,t') = exp [-§ f* dh f* dt 2 (v(ti)v(t2))]. Theory 
becomes less transparent and more cumbersome, but in 
essence remains the same except for the range of very 
strong longitudinal noise. 

Although the longitudinal noise is not effective in 
it cannot be neglected in general. From Eq. i|17|) one 
recognizes that rj z leads to a diffusive precession of the 
transverse components of g. Therefore, the expectation 
values (g±) decay exponentially on the dephasing time 



T <$> 



2 



(21) 



with F z o = dt (rj z (t)r] z (0)) . For short noise cor- 
relations, T0 ~ l/[r n (r] z }]. Thus, longitudinal dephas- 
ing can be neglected during the accumulation period 



\t\ < T"acc = 1/^Tn of transverse noise for (77^) <C v. On 
the other hand, on sufficiently large time scales \t\ 
dephasing always sets in. If the dynamics of the system 
is followed over that long times, the averages of g x and 
g y vanish because of the random precession around the z 
axis. Nevertheless, in this case the transverse amplitude 



g± = (gl + 9 2 y ) 1/2 = (2g+.9-) 1/2 = (2g+g-) 



1/2 



(22) 



which is identical to the absolute value of the off-diagonal 
element of the density matrix, can have a finite expecta- 
tion value. 



III. NOISY DYNAMICS 

In this section, we temporarily neglect regular transi- 
tions (we set b x = 0) and focus on transitions solely due 
to noise. To solve the equations of motion for the propa- 
gation from an initial time to to t > to, we integrate the 
equations of motion Ijl8|l to 



g±(t) = g±(to)±i / dt'fj±(t')g z (t') (23a) 
9z(t) = g z (t )- / dh / dt 2 w(t 1 ,h)g z (h) 

Jto Jto 

dt'[ifj + {t')g-(t ) + cc] (23b) 



with the vertex pair function 

w(h,h) = fj + (h)fj-{t 2 ) + cc 



(24) 



Wherever "cc" occurs, it stands for the complex conju- 
gation of the preceding term. 

Equation (|23bjl may be fed back iteratively into itself 
and into Eq. H23al) . to express g(t) entirely in terms of 
g(to)- Each term arising from this iteration can be visu- 
alized by a diagram. Defining the integral series 



W* (t') = S(t -t')- f dh [ 1 dt 2 w(t u h)S(t 2 -t') 

ft ptl pt2 /**3 

+ dh dt 2 dh dt 4 w(h,h)w(h,h)S(h-t') - ... (25) 

Jto JtQ Jto Jto 



one obtains 

rt' 



g z {t) = J dt'Wt {t') ^g z {h) ~ J t dH*MOM*o)+cc]j . ( 26a ) 
9±(t) = 9±(t ) ± i f dt'n ± {t') f°° dt"Wl' o {t") J g z (t ) - f dt"'[ifj + (t'")g_(t ) + cc] 1 . (26b) 

J t G J — 00 I Jt Q J 
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In the diagrammatic representation, every factor W cor- 
responds to an arbitrary number of vertex pairs w. Cor- 
responding to the definition 1)24)1 . two "polarities" have 
to be considered per pair (cf. Fig. 2Ji) . During the time 
evolution i|2T)|) of g, besides such paired vertices, "excess" 
vertices can appear as first or last vertices during the evo- 
lution period. 



give rise to factors 



(a) 



(b) 



to 



(C ) < + ; n> a a C% i± — 

t t 

FIG. 1: Diagrammatic representation of w and contributions 
to G a .a(t,to). The time axis is chosen to point to the left. 
An encircled + or — represents a noise vertex fj± . (a) A pair 
of vertices connected by a gray bond represents — w{ti,t2). 
A pair of bold circles contains the sum over both possible 
polarities, (b) To leading order in small T n , only diagrams 
with noncrossing vertex contractions contribute. G ZlZ (t, to) 
consists only of noise contraction (wiggly lines) within vertex 
pairs, (c) G+ t +(t,to) consists on a negative excess vertex as 
first vertex (closest to to) and a positive excess vertex as last 
vertex. Contractions can be performed only from pair to pair. 

According to the Wick theorem for Gaussian fields, the 
noise averaging corresponds to all possible contractions 
between noise vertices. The correlations 1)190 imply that 
contractions can be performed only between vertices with 
opposite charge (in a pictorial language, we associate a 
"charge" ±1 with every vertex fj±). Therefore, only neu- 
tral diagrams (where the number of vertices i) + equals 
the number of vertices ?7_) do not vanish. 



A. Dynamics of averages 

Before addressing correlations, we analyze which dia- 
grams contribute to the average (g(t)). Because of the 
neutrality constraint, this average reduces to 



(~g z (t)) = I dt'(Wl{t'))(UU)) 



(27a) 



(g±(t)) = {l- / dt' ^ dt" f dt"> 

L J tQ J — OO J to 

x%(0<V)^(n>}(3 ± (*o)> (27b) 

In factorizing the expectation values we assume that 
noise before and after to is statistically independent. 

To leading order in r n -C T acc , only contractions be- 
tween neighboring vertices (in time order) contribute. In 
(g z (t)), all contractions are performed within the pair 
functions w(t,t') (cf. Fig. QJ)), which (for t — to S> r n ) 



dt'(w(t,t')) = F±(vt) 
with the Fourier transform 

/•OO 

F±(Sl) = / dr cos(Or)/j_(r) . (28) 



Returning to Eq. ()27a)l . we find that the series can be 
easily summed up to 



(§»(*)> = G*,*(t,to)<k(*b)> 



(29a) 



with 



G,,,(t,*o) = / *'«(*')> 

J — OO 

= expj- jf dt'F±(vt')^ . (29b) 

This result was obtained in the Ref. [6j via the differen- 
tial equation 



d(g z (t)) 
dt 



= -F ± (vt)(g z (t)) 



(30) 



which immediately follows from the integral formula 1)29)) 
and confirms the statistical independence of the events 
before and after some moment of time t on a time scale 
much larger than r n . 

Likewise, to leading order in r n <C r acc , (g + (t)) con- 
sists only of contractions linking neighboring pairs due 
to the presence of excess vertices (cf. Fig. QJ;) . Such a 
contraction leads to a factor 



dt'{rj ± (t)fj^(t'))=F±(vt) 



with 



Ff(Q)= / dr e ±tSi V±(r) • 
Jo 

Apparently, F± = F^ + F^ and F^ = F^[* because 
of the even parity of the noise correlator f±(t). Hence, 
resummation of diagrams gives 



with 



(g+(t)) = G+.+ (t,t )(g + {t )) 



pt rOC rt' 

G+,+ (t,t ) = 1- / dt' I dt" I dt'" 

J to J —oo J to 



(31) 



x(v ± (t')Wl(t")^(t'")) 
expj-^ dt'F^(ut')j . (32) 
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In summary, the dynamics of averages can brought into 
the compact form (a = z,±) 

(&»(*)) = Ga,a(t,to){g a (to)} , (33a) 
G a , a (t,to) = e-[*-W-*-(*°)] (33b) 
by means of noise integrals 

0±(t) = / dt'F±(vt') , (34a) 

*/ — OO 

z (t) = f dt'F ± (ut') . (34b) 

•J — OO 

Note that the subscripts of ■& relate to the components of 
the Bloch vector rather than to the noise components. 

B. Dynamics of pair correlations 

Using the diagrammatic approach established above, 
we now proceed to calculate pair correlations of the Bloch 
vector for b x = 0. 



1. Reduction to equal-time correlators 

We start by realizing that 

{g a {t)h(t)) = e-^-^HUmm (35) 

since noise is uncorrelated before and after i, assuming 
(without loss of generality) that t > t. Thus we need to 
calculate only equal-time correlators 

C a[i {t) = (g a (t)g p (t)) . 

The correlations of principal interest are 

(9z(t)g z (t)) = (g z (t)g z (t)) , (36a) 
(9+(t)g-(t)) = (g+(t)g-(t)) • (36b) 

They are invariant under the transformation to the in- 
teraction representation. 

We describe the evolution of the equal-time correlators 

C a p{t) = Gap n s(t, to)C-fs(to) 

in terms of propagators G a p n s(t,to)- In the diagram- 
matic representation, these propagators consist only of 
ladder-like diagrams: The times axes t and t represent 
the legs of the ladder. Noise contractions between the 
sides are the rungs. A contribution of any diagram with 
intersecting or overlapping lines is smaller than that of 
the main sequence by a factor of the order T n /r acc . Since 
only "neutral" diagrams survive averaging, many propa- 
gators vanish. The surviving propagators satisfy "charge 
conservation" 

X (a) + X (0) = x{l) + X(S) 



with the charge function 

{1 for a = + , 

for a = z , 

— 1 for a = — , 

which now refers to propagator indices, not to noise in- 
dices. At the initial time to, the propagator charge is op- 
posite to the charge of the first excess vertex. At the final 
time, the propagator charge coincides with the charge of 
the last excess vertex. 

Among the nonvanishing propagators, several propa- 
gators are mutually dependent. In particular, there is a 
trivial symmetry G a p n s = Gp a ^ corresponding to the 
exchange of the legs of the ladder. Additional relations 
follows from the time-reversal invariance or the complex 
conjugation, such as G z _ jZ _ = G* z+ z+ . 

We therefore find the time evolution of the equal-time 
correlators captured by the relations 

(t, to)C zz {to) 
+2G zz ,+-(t,to)C+-(to) , (37a) 
C + -(t) = G + -^ zz (t,to)C zz (to) 
+[G+_,+_(Mo) 

+G + _,_ + (Mo)]C + _(; ) , (37b) 
C++(t) = G ++ , ++ (t,to)C ++ (to) , (37c) 
C z+ (t) = [G z+ , z+ (t,t ) 

+G z+ , +z (t,to)}C z+ (t ) . (37d) 

It is important to realize that the dynamics of the "neu- 
tral" correlations C ZiZ and C+,_ is decoupled from the 
"non-neutral" C z .± and C± t ±. This implies that the lat- 
ter correlations will not be generated by noise if they are 
absent initially. However, such correlations can be gen- 
erated in general by the transverse components of the 
magnetic field. 

For the convenience of the reader we anticipate the 
results of our subsequent explicit calculation of the prop- 
agators: 



G zz . zz (t, to) 


= j(l + 2e-»). 


(38a) 


Gh — ,zz(t, to) 


-«-*). 


(38b) 


Gzz.H — (t, to) 


= G-| — , zz (t, to) , 


(38c) 


G+-,+-(t, to) 


= \{2 + ie-^+e-^), 


(38d) 


G+--+(t, t ) 


= \{2-ie-^+e-^), 


(38e) 



introducing the abbreviation 

6 a = K(t)-K(to) ■ (39) 

These explicit calculations are based on integral equa- 
tions which we now derive and solve for the propagators 
G ZZtZZ and G zz ^ as two representative examples. The 



diagrammatic derivation of the other propagators is de- 
ferred to the appendix. We refrain from presenting the 
expressions for G++.++, G z +, z +, and G z + t + z since they 
are not relevant for our purposes. 



Gzz,zz(t,to) 



to 



2. Propagator G Z z,zz 

We start with the propagator for the longitudinal cor- 
relator 



G ZZjZZ (£, to) 



oo poo 

" / dt'{W* (t')W* (i')) 

OO J — OO 



dt' 



This expression follows directly from the contribution to 
the autocorrelation of Eq. (|26a[l stemming from the initial 
value g z (t )- 

To leading order in T n /r acc , the contributing diagrams 
have the ladder structure mentioned previously (cf. Fig. 
|2| • On leg segments between two neighboring rungs of 
the ladder, an arbitrary number of noncrossing vertex 
contractions can be performed. Since there are no final 
excess vertices close to time t, all contractions between 
the last rung (at time t\) and the final time t are intra- 
pair contractions. Between the last rung and the last 
but one rung, all contractions are inter-pair contraction, 
forming a ring structure limited by time t\ and t 2 . Then, 
the next contractions on the legs prior to t 2 can be intra- 
pair contractions again. 

Since all possible diagrams prior to t 2 have the same 
structure as the diagrams prior to t, the propagator sat- 
isfies a Bethe-Salpeter integral equation: 

G Z2 , Z2 (Mo) - G 2 z>z (t,t )+2 f dtt f 1 dt 2 Gl z (t,h) 

xf{vti)G +t+ {t u t 2 )G-,-{t u t 2 )f{vt 2 ) 
xG zz , zz (t 2 ,i ) . (40) 

The first term represents the diagram without rungs. In 
the second term, the last pair of rungs at times t\ and 
t 2 is separated from possible additional pairs prior to 
t 2 . On the ring between t\ and £ 2 , the intra-pair noise 
contractions are "oriented". The explicit factor 2 arises 
from the two possible orientations. 

To simplify the integral equation, it is convenient to 
substitute the time variables by variables y = $ z {t) — 
•d z (to) and yi = *& z {ti) ~ $z(to)- This leads to 



G z 



s {y,0) = e~ 2y + 2 p d yi f 1 dy 2 
Jo Jo 

xe- 2 ^-^e-^-^G zz , zz (y 2 ,0) (41) 



This integral equation can now be transformed into a 
differential equation in two steps: 

^-[e 2y G zz , zz (y,0)\ = 2 /* dy 2 e^ 2 

dy Jo 

xG zz , zz (y 2 ,0) (42a) 

-^{e- y ^[e 2y G zz , zz (y,0)]} = 2e»G„,„(y,0)(i42b) 
dy dy 




z (t 2 ,to) 



FIG. 2: Diagrams contributing to G ZZ:ZZ (t,to). First, there 
is the "ladder" without rungs. The diagrams with rungs can 
be represented as propagators G z . z (t,t\) on both legs of the 
ladder, following the ring contraction between the last rungs 
as times ii and t-z. Prior to £2 additional leg propagators 
and rings may occur. In the Bethe-Salpeter equation, these 
additional contraction sum up to G ZZlZZ (t2,to). 



Performing the derivatives, the last equation simplifies to 

^G zz , zz {y, 0) + Z^-G zz , zz {y, 0) = . 
dy dy 

This differential equation must be solved with the initial 
conditions 

G 2Z , ZZ (0,0) = 1 , (43a) 

-^G zz , zz (y,0)| a=0 + 2G ZZ , ZZ (0,0) = 0. (43b) 
dy 

The first one expresses the trivial fact that there is no 
time evolution over a vanishing time interval and follows 
from equation l|4(J|l . The second one stems directly from 
Eq. (|42a|l . Both determine a surprisingly simple form of 
the final result l|38a|l for G zz . zz (t,to). 

Considering an initial state with g(£o) = 9z(t-o) e zy the 
conservation of g 2 implies the normalization 

G zz ,zz + 2G+_ ZZ = 1 . (44) 

This implies the result (|38b|l for G^ , Z z(t,to) which is 

derived also diagrammatically in the appendix. 



3. Propagator G ZZj +_ 

As an example for the diagrammatic calculation of the 
other propagators, we consider here 

/OO ft' />OC ft' 

dt' / dt" / dt' I dt" 
-00 J to J — oo J to 

x(w! a (t')Wl(F)^(t")r ]+ (F')) (45) 
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This expression follows directly from the contribution to 
the autocorrelation of Eq. (|26a[) stemming from the ini- 
tial value g±(to). Its diagrams (cf. Fig- EJ have no excess 
vertices close to the final time t but opposite excess ver- 
tices close to the initial time to- Therefore, close to to 
only inter-pair contractions are possible which are con- 
nected through the first rung at the time t\ . All possible 
contractions between t\ and t have the same structure as 
the diagrams of G ZZjZZ . Therefore, the diagrams can be 
summed implicitly to 



G 



zz,-\-- 



-(*,*o) 



xG +)+ (ti s *o)G-,-(*i,«d) ■ (46) 



After a transition to variables y the integral can be per- 
formed explicitly and yields G zz ^ = G^ ;ZZ . 



G zz , + .(t,t ) 



Gzz,zz(t,ti) 



ti t 



FIG. 3: The propagator G zz ,+- consists of diagrams in which 
the first vertices are excess vertices. They can be contracted 
with a certain number of pairs until the first rung occurs at 
time t\. The contraction of all later pairs yields G zz , zz (£, ti). 

The normalization condition 

G zz , + _(t,t ) + G + _, + _(t,t ) + G + _,_ + (t,t ) = 1 , 

which expresses the conservation of g 2 for an initial trans- 
verse state, determines already the sum of the prop- 
agators G^ (t, to) and G + _, \-(t,to) which enter 

Eqs. i|37|) only in this combination. Nevertheless, both 
propagators G-| ^ (t, to) and G_| . \-(t,to) are calcu- 
lated separately in the appendix. 



IV. TRANSITIONS PRODUCED BY NOISE 
ONLY 

On the basis of the propagators derived in the previ- 
ous section we now evaluate transition probabilities and 
their fluctuations, still focusing on the subcase b x = 0. 
We consider an arbitrary initial state at to = — °o with 
arbitrary Bloch vector g(— oo). We are interested in the 
final time t — oo, where the average is given by 



(fc M) 
(.9±M) 



- e/2 fl±(-oo) 



(47a) 
(47b) 



according to Eqs. dS3>- Since the noise correlators are 
even functions of time differences, only the real parts of 
■d a contribute to the difference t? Q (oo) — d a (— oo). We 
have introduced 

e a = 0,(00) = -<»£ + t£) 

for abbreviation, ft is important to realize that the de- 
cay of the averaged Bloch vector depends only on the 
instantaneous noise correlation. 

Although the final slow amplitudes g± (oo) will be finite 
if the initial amplitudes g±(— oo) were finite, the original 
amplitudes vanish, 



G7±(oo)> =0 



(48) 



due to the diffusive random precession leading to a decay 
on the dephasing time given in Eq. (|21|) . However, 
the transverse amplitude g± decays like g±, 



(0x(oo)) 



- p-e/2 



(49) 



For the variances we obtain - suppressing the obvious 
time arguments (oo, — oo) of the propagators - 



(s 2 M> 

(ffi(oo)) 



— G ZZi 

= \{ S 

= 2G+ 



sflf (-oo) + G zz ,+_gi(-oo) 



,-38 



2^(_oo) - fli(-oo)]}(50a) 



,zz9 z (-oo) 
[G + _ )+ _ + G + _,_+] 5 i(-oo) 



= 3 {2g2 - ' 
and for the fluctuations 



-36 



2^-co)-.g 2 (-oo)]$,Qb) 



([5 5z (oo)] 2 ) = ( ff 2 (co)) - < 3z (co)) : 



^{g 2 + [2e" 3e 

" - 3 V(— )} 



3e 



-20i 2 



&(-oo) 



(51a) 



([^(co)] 2 ) = ( 5 2 (co))-(. 9± (co)) 2 

= i{2g 2 -2e- 30 ff 2 (-co) 



-30 



3e-«] ff i(-co)} . (51b) 



For systems prepared initially in a state with a well de- 
fined g, the correlations C z ± and C±± do not necessarily 
vanish. They are not explicitly considered here since they 
do not couple to G zz and G_| according to Eqs.|23 Fur- 
thermore, an long time scales, phase diffusion due to r\ z 
annihilates these non-neutral correlations anyway. 

To discuss the statistics of transitions between the two 
levels, we consider the system to be located initially in 
the state |f), i.e. g(— oo) = e z . A measurement which 
detects at time t whether the system is in the state |f) 
or | J.) yields the probabilities 



n TT/n (t) = Tri(i ±a z )p(t) 



1 



[l±g z (t)] .(52) 



8 



The trace accounts for the average over quantum fluc- 
tuations for a given realization of noise. The additional 
averaging over the noise leads to the final expectation 
values 

(n TT M> = \[l + e-% (n u (oo)) = \[l-e-®] (53a) 

Since the probabilities n-p-j- /-pj. depend on the noise, they 
are fluctuating quantities themselves. Their fluctuations 
are 

([<ffl u M] 2 ) = 1(1 - 3e- 2e + 2e- 3e ) . (53b) 

Remarkably, the fluctuations of probabilities in gen- 
eral have the same order of magnitude as average val- 
ues. They are weak only if the noise is weak, i.e., 
([(Hl T p(oo)] 2 ) 9 2 /4 for 8 < 1. In the opposite lim- 
iting case 8 3> 1 of strong noise the two levels become 
equipopulated, (IIpp(oo)) = (Iljp(oo)) = |, and the fluc- 
tuation is equal to ([(5I±|p(oo)] 2 ) = 

V. LZ SYSTEM WITH NOISE 

In this section we assume that not only the noise, 
but also the regular part of the Hamiltonian has a finite 
nondiagonal matrix element A corresponding to a finite 
transverse component of the regular magnetic field. As 
it was shown in Ref. 6], if the noise is fast, there exists 
a well-defined time separation for the transition due to 
the noise and due to the regular part of the Hamiltonian. 
Consider matching at times ±i x with £lz C ( x < T"acc- 
The time evolution may be decomposed into three in- 
tervals, (i) From t = — oo to t = —t x one may neglect 
A. The transition is only under the influence of noise. 
Since t x <C r acc , this time evolution is approximately the 
same as from t = — oo to t = in the absence of A. (ii) 
From t = —t x to t = t x the LZ transition occurs. Since 
tx ^> tlz, time evolution is approximately the same as 
from t = — oo to t = oo in the absence of noise, (iii) From 
t = tx to £ = oo the situation is again analogous to the 
regime (i). 

We represent the time evolution operator as 

U{t,t ) = U (t)U(t,t Q )Ul(t Q ) (54) 
and approximate 

U(oo,— oo) = UmUiiUi , (55a) 

Ui = [7(0, -oo) for A = , (55b) 

Uii = £7(oo, — oo) for r] = , (55c) 

U m = E/(oo, 0) for A = (55d) 

for the time evolution operator acting on the slow vector. 

We assume that at t = — oo the density matrix is di- 
agonal (complete decoherence) , i.e., 

g(-oo) = g z (-oo)e z . 



Then 

(g z (-oo)g z (-oo)) = g 2 z (-oo) 

is the only non- vanishing correlator at t = — oo. Since 
T n *C t x , the noise in intervals (i) and (iii) is statisti- 
cally independent and the averaging can be performed 
separately for both intervals. 

A. Interval (i) 

Neglecting A in the first time interval, we find at t = 
—t x the averages 

(&(-*x)> = e- e / 2 . 9z (-oo) , (56a) 
(~9±(-tx)) = , (56b) 

and correlations 

C zz {-tx) = G zz , zz (0, -oo)# 2 (-oo) 

= |(l + 2 e - 3e / 2 )^(_ oo) , (57a) 

C+-{-t x ) = G+_, zz (0,-oo) 5 2 (-oo) 

= i(l-e- 3e / 2 ) 5 2 (-oo). (57b) 

According to Eqs. (|37|l all other independent pair corre- 
lations vanish at t = —t x - 

B. Interval (ii) 

Time evolution at intermediate times from —t x to t x 
is given by 

U& = Ulz ■ 

Hence the density matrix evolves according to 

p(*x) = u hZ p(-t x )ul z 

which implies the transformation 

g(t x ) ■ & = U LZ g(-t x ) ■ &u[ z 

for the transformed Bloch vector g. The explicit trans- 
formation of its components reads 

9+(tx)\ / g+(-t x )\ 

9z(t x ) = U LZ • g z (-t x ) (58a) 
9-(tx)J \9-{-tx)J 

with the rotation matrix 

/ a* 2 -V2a*b* -b* 2 \ 
U LZ = y/2a*b |a| 2 -|6| 2 V2ab* . (58b) 
V -b 2 -V2ab a 2 J 
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From these relations we obtain 

(~g z (tx)) = (fo(-tx)) 

= (H 2 - \b\ 2 ) e - Q / 2 g z (-oo) , (59a) 



(g+(t x )) = U¥ z (§,Hx)) 



-V2a*&*e- e/2 5 ,(-oo) 



(59b) 



for the averages. Analogously, correlations are trans- 
formed according to 



7<5 



(60) 



Since we are ultimately interested only in the correlations 

C zz (po) and C_| (oo) and since the time evolution in 

interval (iii) again follows Eqs. we need to evaluate 
at time t = t x only the two correlations 



C zz (t x ) = (U^) 2 C zz (-t x ) 

LZ 

z + 



+2U^Ut z C+-(-t x ) 



= i[l + 2(1 - 6M 2 H 2 ) e - 30 / 2 ] 5 2 (-^(j6aa) 
C + _(t x ) = U™U™C M (-t x ) 

+[u?ft.u™ + u^Lu^l]c+-(-t x ) 

= |[1-(1- 6|a| 2 |6| 2 ) e - 30 / 2 ] 5 2 (-c^)(61b) 

As a check of these results, we verify the normalization 
relation C zz (t x ) + 2C+_(t x ) = g 2 (-oo). 



C. Interval (iii) 

As already mentioned before, the time evolution in in- 
terval (iii) is absolutely analogous to the evolution in in- 
terval (i). We therefore obtain 

(&(oo)) - G z , z (oo,0){g z (t x )) 

= (|a| 2 H6| 2 )e- e 3z (-oo) , (62a) 
(g+(oo)) = G+ 1+ (oo,Q)(&(* x )) 

xe- e/2 . 9z (-oo) . (62b) 

These results were obtained already in Ref. Q ■ In terms 
of g they mean 

(5,(00)) = (H 2 - \b\ 2 )e- e g z (-^) , (63a) 
(ffx(oo)) = V2|a||6|e- 3e / 4 52 (-oo) . (63b) 

We recall that the averages (g x (oo)) = (g y (oo)) = van- 
ish due to the presence of the longitudinal noise compo- 
nent. The correlations follow from Eqs. H37J) with t = oo 
and to — t x '■ 

C zz (oc) = i[l + 2(l-6|a| 2 |6| 2 )e- 30 ] 5 2 (- & oIfi4a) 
C + -(oo) = i[l-(l-6|a| 2 |6| 2 )e- 30 ] 5 2 (-ao)(.64b) 



In the presence of regular transitions, the transition prob- 
abilities and their fluctuations - given in Eqs. I|53[) in the 
absence of A - finally become 



( n TT/U (°°)) = ^[l±(\a\ 2 -\b\ 2 )e- G ] (65a) 



and 



([OT n/n (oo)] 2 ) 



2|I,|2\„ — 20 
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[l-3(l-4| 0r >|> 



2|i,l2\„— 301 



-2(l-6|a| 2 |6| 2 )e 



(65b) 



To obtain the explicit dependence on the LZ parameter 
one may use \a\ 2 = e~ 2lT7 and the unitarity relation 
|6| 2 = l-|a| 2 . 



VI. RESPONSE SPECTRUM AND INTENSITY 

Let a short pulse of an effective field directed along 
axis i act on the LZ system at some moment of time t. 
Its response at a later moment of t' = t + r is determined 
by the correlation function, Ki t j(t, r) = (si(t).Sj(t + t)). 
It is reasonable to measure the spectral content of the 
response at a fixed time of initial pulse t. It is given by 
the Fourier-transform 



K ia (t,r)e lulT dT 



(66) 



According to general relations l|35l) and l|33l) , in the limit 
t — ► oo such a function vanishes or saturates to a finite 
value. For example, in the case i — j = z and t 3> t x it 
is equal to C zz (t)exp (— J t °° F±(vt')cIt') . Therefore, the 
spectral density of the response contains a 1/w compo- 
nent at small frequency with an intensity that depends 
on the time of excitation. This is not surprising since 
the LZ process violates the homogeneity of time. At infi- 
nite time t this intensity reaches a finite limit C zz (oo). It 
happens since the noise becomes ineffective at sufficiently 
large time and does not randomize the quantum ampli- 
tudes any more. Therefore, any quantum amplitude or 
a component of the Bloch vector g changes after this 
time in an almost deterministic way. On the other hand, 
the spectral intensity remains finite even at t — > oo since 
the value of g reached at the moment, when the noise 
becomes ineffective, is random. 

Another (incoherent) contribution to the spectral den- 
sity arises from the variation of the correlation function 
at a finite time: 



K„(r,t) - K„(oa,t) = K„(pa,t) 



exp / F±(vT')dT' 



(67) 



To estimate the spectral width of this contribution 
we assume that either time t is large enough or the 
noise is weak. Then the exponent in equation l|67|) can 
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be expanded to the linear term. Its spectral width is 
easily estimated as Au)\ = l/r acc . The spectral in- 
tensity of this line at t ~ i acc can be estimated as 
I\ cc C zz (oo) (rfe + r)y) jv. It decreases rapidly with 
growing time t. 

If the LZ processes are repeated randomly with the 
average time r co n between the events (as it happens in 
a gas of colliding atoms or molecules), the 1/w-line is 
smeared out to the width Aujq ~ 1/ r co ii <C A^i , whereas 
its average intensity Iq is proportional to C zz (oo). The 
width of the incoherent line remains the same Awi = 
1/facc, but its intensity is stabilized at a value /i(r co n) oc 
-^^(Tcoii/Tacc) the smaller the larger is r co n [here <p(x) is 
a function describing the decay of the noise correlations, 
which depends on details]. 

Thus, the main contribution to the spectral width of 
the population or induced field fluctuations consists of 
two narrow lines. The first has the width determined by 
collisions and a permanent intensity per particle. The 
second has a permanent width, but its intensity is de- 
termined by collision time. In the rarefied gas, or more 
generally if t co \\ S> r acc , the first line is narrower and 
stronger than the second one. 
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APPENDIX A: DETAILS FOR CORRELATORS 

Here we provide some details on the calculation of cor- 
relation propagators. 



1. Propagator G+_, zz 

The calculation of the propagator G_| %zz can be per- 
formed largely in parallel to the calculation of G zz ^ in 

Sec. IIII B 31 Starting point is the equation 

/>t ft poo poo 

<?+-,„(*,«<>) = dt ' dl ' \ / dt " 

J t() J to J —CO J —CO 

x{fj.(t')fj + (r)wf (t'')w! (¥^ . 

Using the Wick theorem, we perform the contractions 



VII. CONCLUSIONS 

We have calculated the correlation functions of the 
Bloch vector components in the LZ process subject to 
a fast Gaussian noise. These correlators determine the 
linear response of the LZ system to a weak, time depen- 
dent probe signal. Two-time correlation functions are 
factorizable at the time scale much larger than the noise 
correlation time r n due to statistical independence of ran- 
dom processes. Thus the main problem is the calculation 
of simultaneous averages. The condition of the fast noise 
allows to use for these averages the ladder graphs only. 
The resulting Bethe-Salpeter equation can be reduced to 
differential equations which are exactly solvable. In gen- 
eral, fluctuations are strong, of the same order of magni- 
tude as the average values. 

Since the noise is ineffective for transitions at suffi- 
ciently large time t 3> r acc , the transition probabilities or 
the components of the Bloch vector remain determinis- 
tically coherent after this time. Therefore, the spectrum 
of the fluctuations contains a narrow 1 /u> line, whose in- 
tensity is determined by the distribution of these values 
at the time when the noise became ineffective. Besides 
of this line, there exists another narrow line with a width 
of about 1 / r acc and an intensity depending on the pulse 
time and going to zero when this time goes to infinity. 
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G+_,„(t,td) = / &iG + , + (t,ti)G_,_(t,ti) 
y.f{vt 1 )G zz . zz {t ll t ) , 

where t\ is the time of the last rung of the ladder (last 
means closest to t, see Fig.0J. After a substitution from 
t\ to y\ the integral is elementary with the result 

G + _,„(t,t ) = i(l-e- 3 ^W-^( t »)]) . 



<— 1 

G + ., zz (t,t ) 

M-^ 


iz 


t 


to 




:Z 




z(tl,to) 



FIG. 4: The diagrams of , zz are analogous to the dia- 
grams of G zz ,+- with the only difference that now the excess 
vertices are the last vertices. 
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2. Propagator Gh ,-| 

From the equations of motion (|26|l one extracts 
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pt pOO pt 

G + _,+_(Mo) = 1- dt' dt" / dt"'(^(t')W*' (t")fj + (t'")) 

J to J — oo ./in 



'to 

- [dp 

ft 

dt' 

'to 



to 

t" 



rff / d?"(fj + (?)Wl(F')fi-(i">)) 

J t 

pt" pt poo pi" 

dt" / dt'" / dt' / dt" / 

) J to J to J — oo ./ to 



After averaging, the sum over the ladder diagrams can be captured by a Bethe-Salpeter equation 
G+_,+_(Mo) - G + , + (Mo)G_,_(Mo) 

+ / dti / dt 2 G+, + (*,*i)G_ 1 _(t,ti)/(i^i)G„,„(t 1 ,t a )/(i^a)G +1+ (*2,*o)G_ 1 _(<a,to) 

Jt Jt 

where ii is the time of the last rung and £2 the time of the first rung (see Fig. EJ. All contractions in between sum 
up to G zz , zz (ti 7 i 2 ). Since the latter propagator is explicitly known, the integrals can be performed directly, yielding 
Eq. i38dl) ' 



<i 

G + 




, + -(t,to) 


.+ 


t 




to 


< + ; n <c\ 




.+ 


J £\ £1 






t 


+ 


to 


e-)i 


it—e <jf™\ 


.+ 



G zz , zz (t 1; t 2 ) 
-it »-i e-e- 



to 



FIG. 5: G-| ,_| consists of diagrams with excess vertices as first and last vertices on both time legs. There are disconnected 

diagrams (where the "ladder" has no rungs) and connected diagrams (where the ladder has rungs). In the second case, the excess 
vertices are are connected by rungs as ti and t^ and between t\ and £2 additional pairs can be contracted as in Gz Z ,zz(ti,t2). 



3. Propagator , ^ 

Similar to the previous propagator, one now has 

pt />oo pt pt poo pt 

G + _,_ + (i,fo) = dt' dt" / dt'" / dt' / dt" / (fi"'<<J_(f)<'(t /, )<7-(* / ")^+(«')<(«")^-(*" , )> ■ 

J to J — OO J to J to J —OO J to 

The configuration of the excess vertices necessitates at least two rungs at times t\ and £2 (see. Fig.[SJ). The contractions 
in between again sum up to G ZZtZZ (ti, i2), leading to 

G + _,_+(Mo) = f f tft 2 G + , + (t,ti)(3_,_(t,ti)/(^i)G„,„(ti,t2)/(^2)G+,+ (t2,to)G_,_(t2,t^ 

J tQ J to 
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G + .,. + (t,t ) 


+ 




t 

,0, . 




1 

... /-\ 







-a a; 


Mi 

G zz ,zz(tl,t 2 ) 
$-n 




-A . 


+ 



t tl t 2 t 



FIG. 6: The diagrams of , \. as similar to these of 

G+_,+- with the difference that the polarity of the first ex- 
cess vertices is inverted. Therefore, no disconnected diagrams 
can be formed. 



and eventually to Eq. (|38el) . 
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